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Abstract 

A simple scheme for communication over MIMO broadcast channels is introduced which 
adopts the lattice reduction technique to improve the naive channel inversion method. Lattice basis 
reduction helps us to reduce the average transmitted energy by modifying the region which includes 
the constellation points. Simulation results show that the proposed scheme performs well, and as 
compared to the more complex methods (such as the perturbation method [1]) has a negligible 
loss. Moreover, the proposed method is extended to the case of different rates for different users. 
The asymptotic behavior (SNR — > oo) of the symbol error rate of the proposed method and the 
perturbation technique, and also the outage probability for the case of fixed-rate users is analyzed. It 
is shown that the proposed method, based on LLL lattice reduction, achieves the optimum asymptotic 
slope of symbol-error-rate (called the precoding diversity). Also, the outage probability for the case 
of fixed sum-rate is analyzed. 



I. Introduction 

In the recent years, communications over multiple- antenna fading channels has attracted 
the attention of many researchers. Initially, the main interest has been on the point-to-point 
Multiple-Input Multiple-Output (MIMO) communications [2]-[6]. In [2] and [3], the authors 
have shown that the capacity of a MIMO point-to-point channel increases linearly with the 
minimum number of the transmit and the receive antennas. 

More recently, new information theoretic results [7], [8], [9], [10] have shown that 
in multiuser MIMO systems, one can exploit most of the advantages of multiple-antenna 
systems. It has been shown that in a MIMO broadcast system, the sum-capacity grows 
linearly with the minimum number of the transmit and receive antennas [8], [9], [10]. To 
achieve the sum capacity, some information theoretic schemes, based on dirty-paper coding, 
are introduced. Dirty-paper coding was originally proposed for the Gaussian interference 
channel when the interfering signal is known at the transmitter [11]. Some methods, such 
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as using nested lattices, are introduced as practical techniques to achieve the sum-capacity 
promised by the dirty -paper coding [12]. However, these methods are not easy to implement. 

As a simple precoding scheme for MIMO broadcast systems, the channel inversion 
technique (or zero-forcing beamforming [7]) can be used at the transmitter to separate the data 
for different users. To improve the performance of the channel inversion technique, a zero- 
forcing approximation of the dirty paper coding (based on QR decomposition) is introduced 
in [7] (which can be seen as a scalar approximation of [12]). However, both of these methods 
are vulnerable to the poor channel conditions, due to the occasional near-singularity of the 
channel matrix (when the channel matrix has at least one small eigenvalue). This drawback 
results in a poor performance in terms of the symbol-error-rate for the mentioned methods 
[1]. 

In [1], the authors have introduced a vector perturbation technique which has a good 
performance in terms of symbol error rate. Nonetheless, this technique requires a lattice 
decoder which is an NP-hard problem. To reduce the complexity of the lattice decoder, in 
[13]— [16], the authors have used lattice-basis reduction to approximate the closest lattice 
point (using Babai approximation). 

In this paper, we present a transmission technique for the MIMO broadcast channel 
based on the lattice-basis reduction. Instead of approximating the closest lattice point in the 
perturbation problem, we use the lattice-basis reduction to reduce the average transmitted 
energy by reducing the second moment of the fundamental region generated by the lattice 
basis. This viewpoint helps us to: (i) achieve a better performance as compared to [14], (ii) 
expand the idea for the case of unequal-rate transmission, and (iii) obtain some analytic results 
for the asymptotic behavior (SNR — > oo) of the symbol-error-rate for both the proposed 
technique and the perturbation technique of [1]. 

The rest of the paper is organized as the following: Sections II and III briefly describe 
the system model and introduce the concept of lattice basis reduction. In section IV, the 
proposed method is described and in section V, the proposed approach is extended for the 
case of unequal-rate transmission. In section VI, we consider the asymptotic performance of 
the proposed method for high SNR values, in terms of the probability of error. We define the 
precoding diversity and the outage probability for the case of fixed-rate users. It is shown 
that by using lattice basis reduction, we can achieve the maximum precoding diversity. For 
the proof, we use a bound on the orthogonal deficiency of an LLL-reduced basis. Also, an 
upper bound is given for the probability that the length of the shortest vector of a lattice 
(generated by complex Gaussian vectors) is smaller than a given value. Using this result, 
we also show that the perturbation technique achieves the maximum precoding diversity. 
In section VII, some simulation results are presented. These results show that the proposed 
method offers almost the same performance as [1] with a much smaller complexity. As 
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compared to [14], the proposed method offers almost the same performance. However, by 
sending a very small amount of side information (a few bits for one fading block), the 
modified proposed method offers a better performance with a similar complexity. Finally, in 
section VIII, some concluding remarks are presented. 

II. System Model and Problem Formulation 

We consider a multiple-antenna broadcast system with N t transmit antennas and N r 
single- antenna users (N t > N r ). Consider y = [y t , ...,y Nr ] T , x = [x\, x Nt ] T , w = 
[wi, wn t ] t , and the iV r x N t matrix H, respectively, as the received signal, the transmitted 
signal, the noise vector, and the channel matrix. The transmission over the channel can be 
formulated as, 

y = Hx + w. (1) 

The channel is assumed to be Raleigh, i.e. the elements of H are i.i.d. with the zero- 
mean unit-variance complex Gaussian distribution and the noise is i.i.d. additive Gaussian. 
Moreover, we have the energy constraint on the transmitted signal, E(||x|| 2 ) = 1. The energy 
of the additive noise is a 2 per antenna, i.e. E(||w|| 2 ) = N r a 2 . The Signal-to-Noise Ratio 
(SNR) is defined as p = 4y. 

In a broadcast system, the receivers do not cooperate with each other (they should 
decode their respective data, independently). The main strategy in dealing with this restriction 
is to apply an appropriate precoding scheme at the transmitter. The simplest method in this 
category is using the channel inversion technique at the transmitter to separate the data for 
different users: 

s = H+u, (2) 
where H + = H H (HH H )~ 1 , and H H is the Hermitian of H. Moreover, s is the transmitted 

g 

signal before the normalization (x = — , is the normalized transmitted signal), and 

u is the data vector, i.e. Ui is the data for the i'th user. For N t = N r (the number of transmit 
antennas and the number of users are equal), the transmitted signal is 

s = H^u. (3) 

The problem arises when H is poorly conditioned and ||s|| becomes very large, resulting 
in a high power consumption. This situation occurs when at least one of the singular values of 
H is very small which results in vectors with large norms as the columns of H + . Fortunately, 
most of the time (especially for high SNRs), we can combat the effect of a small singular 
value by changing the supporting region of the constellation which is the main motivation 
behind the current article. 
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When the data of different users are selected from Z[z], the overall constellation can 
be seen as a set of lattice points. In this case, lattice algorithms can be used to modify 
the constellation. Especially, lattice-basis reduction is a natural solution for modifying the 
supporting region of the constellation. 

III. Lattice-Basis Reduction 

Lattice structures have been frequently used in different communication applications 
such as quantization or decoding of MIMO systems. A real (or complex) lattice A is a 
discrete set of A-D vectors in the real Euclidean space M. N (or the complex Euclidean space 
C^) that forms a group under ordinary vector addition. Every lattice A is generated by the 
integer linear combinations of some set of linearly independent vectors bi, ■ • • , \>m G A, 
where the integer M, M < A, is called the dimension of the lattice A. The set of vectors 
{bi, ■ ■ • , h M } is called a basis of A, and the matrix B = [b 1; ■ • • , b M ], which has the basis 
vectors as its columns, is called the basis matrix (or generator matrix) of A. 

The basis for representing a lattice is not unique. Usually a basis consisting of relatively 
short and nearly orthogonal vectors is desirable. The procedure of finding such a basis for 
a lattice is called Lattice Basis Reduction. A popular criterion for lattice-basis reduction is 
to find a basis such that ||bi|| ■ ... ■ ||bji/|| is minimized. Because the volume of the lattice§ 
does not change with the change of basis, this problem is equivalent to minimizing the 
orthogonality defect which is defined as 

A (iib 1 |i 2 iib 2 |i 2 ...||b M in 

6 = — ™ — • (4) 

The problem of finding such a basis is NP-hard [17]. Several distinct sub-optimal 
reductions have been studied in the literature, including those associated to the names 
Minkowski, Korkin-Zolotarev, and more recently Lenstra-Lenstra and Lovasz (LLL) [18]. 

An ordered basis (bi, ■ ■ • , b M ) is a Minkow ski-Reduced Basis [19] if 

• hi is the shortest nonzero vector in the lattice A, and 

• For each k = 2, M, is the shortest nonzero vector in A such that (bi, ■ • ■ , b^) 
may be extended to a basis of A. 

Minkowski reduction can be seen as a greedy solution for the lattice-basis reduction problem. 
However, finding Minkowski reduced basis is equivalent to finding the shortest vector in the 
lattice and this problem by itself is NP-hard. Thus, there is no polynomial time algorithm 
for this reduction method. 

In [20], a reduction algorithm (called LLL algorithm) is introduced which uses the 
Gram-Schmidt orthogonalization and has a polynomial complexity and guarantees a bounded 

2 Volume of the lattice generated by B is (det B H B) 5 . 
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orthogonality defect. For any ordered basis of A, say (bi, • ■ ■ , l>Af)» one can compute an 
ordered set of Gram-Schmidt vectors, Nc>i, ■ • • , b A /j, which are mutually orthogonal, using 
the following recursion: 

bi = hi - Y^iVijbj, with 

< b t ,bj > (5) 

< b i> b i > 

where < •, ■ > is the inner product. An ordered basis (b 1; ■ • • , h M ) is an LLL Reduced Basis 
[20] if, 

• \\pi>ij\\ < | for 1 < i < j < M, and 

• V l|bi|| 2 < ||b i+ i + Hi +lti bi\\ 2 

where | < p < 1, and ^b 1; ■ ■ • , b M J is the Gram-Schmidt orthogonalization of the ordered 

basis and bj = Yl)=i f^ij^j f° r & = 1, Af. 

It is shown that LLL basis-reduction algorithm produces relatively short basis vectors 
with a polynomial-time computational complexity [20]. The LLL basis reduction has found 
extended applications in several contexts due to its polynomial-time complexity. In [21], 
LLL algorithm is generalized for Euclidean rings (including the ring of complex integers). 
In this paper, we will use the following important property of the complex LLL reduction 
(for p = §): 

Theorem 1 (see [21]): Let A be an M-dimensional complex lattice and B = [bi...b&f] 
be the LLL reduced basis of A. If 5 is the orthogonality defect of B, then, 

V6<2 M W- 1 \ (6) 

IV. Proposed Approach 

Assume that the data for different users, Ui, is selected from the points of the integer 
lattice (or from the half-integer grid [22]). The data vector u is a point in the Cartesian 
product of these sub-constellations. As a result, the overall receive constellation consists 
of the points from I? Nr , bounded within a 2A r r -dimensional hypercube. At the transmitter 
side, when we use the channel inversion technique, the transmitted signal is a point inside 
a parallelotope whose edges are parallel to vectors, defined by the columns of H + . If the 
data is a point from the integer lattice I? Nr , the transmitted signal is a point in the lattice 
generated by H + . When the squared norm of at least one of the columns of H + is too large, 
some of the constellation points require high energy for the transmission. We try to reduce 
the average transmitted energy, by replacing these points with some other points with smaller 
square norms. However, the lack of cooperation among the users imposes the restriction that 
the received signals should belong to the integer lattice Z 27Vr (to avoid the interference among 
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the users). The core of the idea in this paper is based on using an appropriate supporting 
region for the transmitted signal set to minimize the average energy, without changing the 
underlying lattice. This is achieved through the lattice-basis reduction. 

When we use the continuous approximation (which is appropriate for large constella- 
tions), the average energy of the transmitted signal is approximated by the second moment 
of the transmitted region [22]. When we assume equal rates for the users, e.g. R bits per 
user (-j bits per dimension), the signal points (at the receiver) are inside a hypercube with 
an edge of length a where 

a = 2 R ' 2 . (7) 

Therefore, the supporting region of the transmitted signal is the scaled version of the fun- 
damental region of the lattice generated by H + (corresponding to its basis) with the scaling 
factor a. Note that by changing the basis for this lattice, we can change the corresponding 
fundamental region (a parallelotope generated by the basis of the lattice and centered at 
the origin). The second moment of the resulting region is proportional to the sum of the 
squared norms of the basis vectors (see Appendix A). Therefore, we should try to find a 
basis reduction method which minimizes the sum of the squared norms of the basis vectors. 
Figure 1 shows the application of the lattice basis reduction in reducing the average energy 
by replacing the old basis with a new basis which has shorter vectors. In this figure, by 
changing the basis ai,a 2 (columns of H + ) to t>i,b 2 (the reduced basis), the fundamental 
region T, generated by the original basis, is replaced by J 7 ', generated by the reduced basis. 

Among the known reduction algorithms, the Minkowski reduction can be considered 
as an appropriate greedy algorithm for our problem. Indeed, the Minkowski algorithm is the 
successively optimum solution because in each step, it finds the shortest vector. However, 
the complexity of the Minkowski reduction is equal to the complexity of the shortest-lattice- 
vector problem which is known to be NP-hard [23]. Therefore, we use the LLL reduction 
algorithm which is a suboptimum solution with a polynomial complexity. 

Assume that B = H + U is the LLL-reduced basis for the lattice obtained by H + , 
where U is an N r x N r unimodular matrix (both U and U _1 have integer entries). We use 
x = Bu' = H + Uu' as the transmitted signal where 

u' = U _1 u mod a (8) 

is the precoded data vector, u is the original data vector, and a is the length of the edges 
of the hypercube, defined by ©. At the receiver side, we use modulo operation to find the 
original data: 

y = Hx + n = HH+U(U x u mod a) + n = I^IT 1 !! mod a) + n (9) 
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Fig. 1. Using lattice-basis reduction for reducing the average energy 



= UU u mod a + n = u mod a + n. (10) 

In obtaining (flOl) from ©, we use the fact that U and U^ 1 have integer entries. 

In this method, at the beginning of each fading block, we reduce the lattice obtained 
by H + and during this block the transmitted signal is computed using <£8j) . Neglecting the 
preprocessing at the beginning of the block (for lattice reduction), the complexity of the 
precoding is in the order of a matrix multiplication and a modulo operation. Therefore, the 
complexity of the proposed precoding method is comparable to the complexity of the channel 
inversion method. However, as we will show by the simulation results, the performance of 
this method is significantly better, and indeed, is near the performance of the perturbation 
method, presented in [1]. 
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In the perturbation technique [1], the idea of changing the support region of the con- 
stellation has been implemented using a different approach. In [1], u' = u + al is used as the 
precoded data, where the integer vector 1 is chosen to minimize ||H + (u+al) || . This problem is 
equivalent to the closest-lattice-point problem for the lattice generated by aH + (i.e. finding 
the lattice point which is closer to — H + u). Therefore, in the perturbation technique, the 
support region of the constellation is a scaled version of the Voronoi region [24] of the 
lattice. In the proposed method, we use a parallelotope (generated by the reduced basis of 
the lattice), instead of the Voronoi region. Although this approximation results in a larger 
second moment (i.e. higher energy consumption), it enables us to use a simple precoding 
technique, instead of solving the closest-lattice-point problem. 

For the lattice constellations, using a parallelotope instead of the Voronoi region (pre- 
sented in this paper) is equivalent with using the Babai approximation instead of the exact 
lattice decoding (previously presented in [14]). The only practical difference between our 
lattice-reduction-aided scheme and the scheme presented in [14] is that we reduce H + , 
while in [14], H H is reduced. As can be observed in the simulation results (figure 2), this 
difference has no significant effect on the performance. However, the new viewpoint helps 
us in extending the proposed method for the case of variable-rate transmission and obtaining 
some analytical results for the asymptotic performance. 

The performance of the proposed lattice-reduction aided scheme can be improved by 
combining it with other schemes, such as regularization [25] or the V-BLAST precoding 
[26], or by sending a very small amount of side information. In the rest of this section, we 
present two of these modifications. 

A. Regularized lattice-reduction-aided precoding 

In [25], the authors have proposed a regularization scheme to reduce the transmitted 
power, by avoiding the near- singularity of H. In this method, instead of using H + = 
H H (HH")" 1 , the transmitted vector is constructed as 

x = H H (HH H + al) _1 u (11) 

where a is a positive number. To combine the regularized scheme with our lattice-reduction- 
aided scheme, we consider B r as the matrix corresponding to the reduced basis of the 
lattice generated by H H (HH H + al) 1 . When we use the regularization, the received signals 
of different users are not orthogonal anymore and the interference acts like extra noise. 
Parameter a should be optimized such that the ratio between the power of received signal 
and the power of the effective noise is minimized [25]. 
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B. Modified lattice-reduction-aided precoding with small side information 

In practical systems, we are interested in using a subset of points with odd coordinates 
from the integer lattice. In these cases, we can improve the performance of the proposed 
method by sending a very small amount of side information. When the data vector u 
consists of odd integers, using the lattice-basis reduction may result in points with some even 
coordinates (i.e. U _1 u has some even elements), instead of points with all-odd coordinates 
in the new basis. For this case, in ([8]), the set of precoded data u' is not centered at the 
origin, hence the transmitted constellation (which includes all the valid points Bu') is not 
centered at the origin. Therefore, we can reduce the transmitted energy and improve the 
performance by shifting the center of the constellation to the origin. It can be shown that 
the translation vector is equal to (U -1 [l + i, 1 + i, • ■ ■ , 1 + i] T + [1 + i, 1 + i, ■ ■ ■ , 1 + i] T ) 
mod 2 where i = a/— 1. When we use this shifted version of the constellation, we must 
send the translation vector to the users (by sending 2 bits per user) at the beginning of the 
block. However, compared to the size of the block of data, the overhead of these two bits is 
negligible. 

The above idea of using a shift vector can be also used to improve the perturbation 
technique (if we only use the odd points of the lattice). After reducing the inverse of the 
channel matrix and obtaining the bits (corresponding to the shift vector) at the beginning of 
each fading block, the closest point to the signal computed in equation ([8]) can be found by 
using the sphere decoder. Then, the transmitted signal is obtained by 

x = B (u + a\ + u par ) , (12) 

where u par is the zero-one shift vector, which is computed for users at the beginning of the 
fading block, and the perturbation vector 1 is an even integer vector such that the vector x 
has the minimum energy. This method which can be considered as modified perturbation 
method outperforms the perturbation method in [1]. When we are not restricted to the odd 
lattice points, using ([121) instead of H + (u' + al) does not change the performance of the 
perturbation method. It only reduces the complexity of the lattice decoder [27]. 

V. Unequal-Rate Transmission 

In the previous section, we had considered the case that the transmission rates for 
different users are equal. In some applications, we are interested in assigning different rates 
to different users. Consider Rx, ...,R Nr as the transmission rates for the users (we consider 
them as even integer numbers). Equation ([8]) should be modified as 

u' = U _1 u mod a, (13) 
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where the entries of a = [at, aN r ] T are equal to 

a, = 2 R ^. (14) 
Also, at the receiver side, instead of Q and (flOl) . we have 

y = Hx + n = HH + U(U x u mod a) + n = l^U^u mod a) + n (15) 
= UU x u mod a+n = u mod a + n. (16) 

If we are interested in sum-rate, instead of individual rates, we can improve the perfor- 
mance of the proposed method by assigning variable rates to different users. We assume that 
the sum-rate (rather than the individual rates) is fixed and we want to reduce the average 
transmitted energy. To simplify the analysis, we use the continuous approximation which has 
a good accuracy for high rates. 

Considering continuous approximation, the sum-rate is proportional to the logarithm of 
the volume of the lattice with basis B and the average energy is proportional to the second 
moment of the corresponding prallelotope, which is proportional to Ylf=i II ^ II 2 = tr BB H 
(see Appendix A). The goal is to minimize the average energy while the sum-rate is fixed. 
We can use another lattice generated by B' with the same volume, where its basis vectors 
are scaled versions of the vectors of the basis B, according to different rates for different 
users. Therefore, we can use B' = BD instead of B (where D is a unit determinant N r x N r 
diagonal matrix which does not change the volume of the lattice). For a given reduced basis 
B, the product of the squared norms of the new basis vectors is constant: 

||b' 1 || 2 ||b' 2 || 2 ...||b' iVr || 2 = (||b 1 || 2 ||b 2 || 2 ...||b Arr || 2 )detD 

= ||b 1 || 2 ||b 2 || 2 ...||b^|| 2 = const. ^ 

The average energy corresponding to the new lattice basis should be minimized. When 
we use the modified basis B' instead of B, the average energy is proportional to Y2i=i II || 2 = 
trB'B /H (see Appendix A). According to the arithmetic-geometric mean inequality, YliZi II II 2 
trB'B /H is minimized iff 

||bi|| = |m|| = ... = ||bi rp ||. (18) 

Therefore, 

mintrB'B /H = N r (||b , 1 || 2 ||b 2 || 2 ...||b / iVr || 2 )^ = N r (||b 1 || 2 ||b 2 || 2 ...||b JV J 2 )^ (19) 

Having the matrix B, the columns of matrix B' can be found using the equation (fT8l) 
and trB'B /H can be obtained by (fT9l) . Now, for the selection of the reduced basis B, we should 
find B such that ||b 1 || 2 ||b 2 || 2 ...||b7v r || 2 is minimized. Because det B H B = det(H + ) H H + is 

lib ll 2 llb II 2 lib II 2 

given, the best basis reduction is the reduction which maximizes J dct b' h b| ' or m 
other words, minimizes the orthogonality defect. 
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In practice, we use discrete values for the rate, and sometimes, we should assign the rate 
zero to some users (when their channel is very bad). In this case, for the rate assignment for 
other users, we use the lattice reduction on the corresponding sublattice. It should be noted 
that the average transmit power is fixed per channel realization and no long-run averaging 
is considered, and no long-run power allocation is used. 

VI. Diversity and Outage Probability 

In this section, we consider the asymptotic behavior (p — > oo) of the symbol error 
rate (SER) for the proposed method and the perturbation technique. We show that for both 
of these methods, the asymptotic slope of the SER curve is equal to the number of transmit 
antennas. By considering the outage probability of a fixed-rate MIMO broadcast system, 
we will show that for the SER curve in high SNR, the slope obtained by the proposed 
method has the largest achievable value. Also, we analyze the asymptotic behavior of the 
outage probability for the case of fixed sum-rate. We show that in this case, the slope of 
the corresponding curve is equal to the product of the number of transmit antennas and the 
number of single-antenna users. 

A. Fixed-rate users 

When we have the Channel-State Information (CSI) at the transmitter, without any 
assumption on the transmission rates, the outage probability is not meaningful. However, 
when we consider given rates Ri,...,Rn t for different users, we can define the outage 
probability P out as the probability that the point R Nr ) is outside the capacity region. 

Theorem 2: For a MIMO broadcast system with N t transmit antennas, N r single- 
antenna receivers (Nt > N r ), and given rates Ri, Rjy r , 

lim Z^lf^l < Nt _ (20) 

p^oo log 

Proof: Define P out \ as the probability that the capacity of the point-to-point system 
corresponding to the first user (consisting of N t transmit antennas and one receive antenna 
with independent channel coefficients and CSI at the transmitter) is less than R^. 

P outl = Pr{log (1 + p\\YL X f) < Rt} (21) 

where hi is the vector defined by the first row of H. Note that the entries of hi have iid 
complex Gaussian distribution with unit variance. Thus, its square norm has a chi square 
distribution. We have, 

Pr{log(2p||hi|| 2 ) <i?i} (22) 
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(N t -1)\ 

We are interested in the large values of p. For p > 2 Rl ~ l , 



( 2 Rl 1 

PrjllM^-j (23) 

2P f llhlll 2(x)dx (24) 

o 

x Nt ~ x t~ x dx (25) 



2 K 1 2 K 1 

a^ 1 -^-* > [ 2 " — 1 ./ '^e" 1 (26) 



e 1 / 2 P 



o (^t-i)i -Jo (Nt-iy. 

x Nt - x dx (27) 

(28) 



(N t - 1)1 Jo 

2 NtRi c 



where c = 2 ^ t ^, is a constant number. Now, 



log (1 + y0 ||h 1 || 2 ) <log(2p||h 1 || 2 ) for P>^rV ^ 

\\"-i\\ 

^ Um -logPrjlogll + pHhJ 2 )^^} 

p^oo log p 

< lim -logPr{log(2p||h 1 || 2 )<i? 1 } (3i) 

p^oo log p 

-log— — 

< lim p -^— = N t (32) 

p^oo log p 

=j. Um ~ bgjP °" fl < iV t . (33) 

p->co log p 

According to the definition of P ou ti, Pout > Pouti- Therefore, 

lim -hsP*** < lim -log P^i < ^ (34) 

p->oo log p P^oo log p 

■ 

We can define the diversity gain of a MIMO broadcast constellation or its precoding 
- log P e 

diversity as lim^oo — where P e is the probability of error. Similar to [28, lemma 5], 

logp 
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we can bound the precoding diversity by lim^oo — - — . Thus, based on theorem 2, the 



log Po ut 

logp 



maximum achievable diversity is N t . 

We show that the proposed method (based on lattice-basis reduction) achieves the 
maximum precoding diversity. To prove this, in lemma 1 and lemma 2, we relate the length 
of the largest vector of the reduced basis B to ^h h (the minimum distance of the lattice 
generated by H H ). In lemma 3, we bound the probability that g?h h is too small. Finally, 
in theorem 3, we prove the main result by relating the minimum distance of the receive 
constellation to the length of the largest vector of the reduced basis B, and combining the 
bounds on the probability that dip is too small, and the probability that the noise vector is 
too large. 

Lemma 1: Consider B = [t>i...b/vf] as an N x M matrix, with the orthogonality defect 
5, and B~ H = [ai.-.a^f] as the inverse of its Hermitian (or its pseudo-inverse if M < N). 
Then, 

max{ bi \\b M } < . — n n ^ (35) 

mm{||ai||,..., ||a M ||} 

and 

max{ ai , a M } < . f|| ,, ,, ,, r (36) 

min{Nbi||, ||b M ||} 

Proof: Consider bj as an arbitrary column of B. The vector bj can be written as 
b. + Ylii^j c i,fov wnere b- is orthogonal to hj for i ^ j. Now, [b 1 ...b i _ 1 b-b i+1 ...bj\/] can be 
written as BP where P is a unit-determinant M x M matrix (a column operation matrix): 

||b 1 || 2 ...||b J _ 1 || 2 .||b J || 2 .||b m || 2 ...||b A/ || 2 (37) 

= 5 det B H B = 5 det P H B H BP (38) 

= 5 det ([b 1 ...b i _ 1 b^b l+1 ...bA/] H [b 1 ...b i „ 1 b^b i+1 ...b A/ ]) . (39) 

According to the Hadamard theorem: 

det ([b 1 ...b i _ 1 b^b l+1 ...b A/ ] H [b 1 ...b J _ 1 b;b i+1 ...b A/ ]) < (40) 

iib 1 || 2 ...||b J _ 1 || 2 .||b:ii 2 .||b l+1 || 2 ...||b A/ || 2 . (4D 

Therefore, 

||b 1 || 2 ...||b J _ 1 || 2 .||b J || 2 .||b m || 2 ...||b A/ || 2 <5||b 1 || 2 ...||b J _ 1 || 2 .||b^|| 2 .||b m || 2 ...||b A/ || 2 (42) 

=> \\hi\\ < V6\\m (43) 
Also, B + B = I results in <aj, bj> = 1 and <a i; hj> = for i ^ j. Therefore, 

1 = <a i; b,> = <a l ,(b; + ^c iij b i )> = <ai,b> (44) 
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Now, aj and b' { , both are orthogonal to the (M — 1) -dimensional subspace generated by the 
vectors bj (j ^ i). Thus, 

1= <a i; b> = Haill-llb-ll > ||ai||.JM- (45) 

=* i > llb.il M m 

=> INI < n4 (47) 

1 1 a i 1 1 

The above relation is valid for every i, 1 < % < M. Without loss of generality, we can 
assume that max-dlb^l, ||b A /||} = ||b fc ||: 

max{||bi||, ||b M ||} = llbjtll < (48) 

||a fc || 

< ^ . (49) 

mm{||ai||, ||a M ||} 

Similarly, by using (1471) . we can also obtain the following inequality: 

V6 

max{ ai a M } < . — (50) 

mm{||bi||,..., ||b M ||} 

■ 

Lemma 2: Consider B = [bi...bj\/] as an LLL-reduced basis for the lattice generated 
by H + and d H H as the minimum distance of the lattice generated by H H . Then, there is a 
constant olm (independent of H) such that 

max{||b 1 ||,...,||b M ||}<^. (51) 
Proof: According to the theorem 1, 

^< 2 M(M-l)_ (52) 

Consider B H = [a 1; ...,a M ]- By using lemma 1 and (l52l) . 

^ 2 m(m-i) 

max{||bi||,..., ||b M ||} < —t-ft, — n n ttt < —r-rr, — n n nrr (53) 

min{||ai||, ||a M ||} mm{||ai||, ||a M ||j 

The basis B can be written as B = H + U for some unimodular matrix U: 

B H = ((H + U) H ) + = (U H BT H ) + = H H tT H . (54) 

Noting that U _H is unimodular, B~ H = [a 1; a M ] is another basis for the lattice generated 
by H H . Therefore, the vectors ai, ...,3lm are vectors from the lattice generated by H H , and 
therefore, the length of each of them is at least dip: 

||aj >d n » for 1 < % < M (55) 
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min{||a 1 1|, ||a w ||} > d n » (56) 

nM(M-l) 

(USD and (|56]) =j> max{||bi||, ||b M ||} < — • (57) 

■ 

Lemma 3: Assume that the entries of the N x M matrix H has independent complex 
Gaussian distribution with zero mean and unit variance and consider <i H as the minimum 
distance of the lattice generated by H. Then, there is a constant (3n,m such that 

(3 N , M s 2N for M<N 

Pr {dn ^ 2JV ^ { _ (ln£) AT +1) 1} f or M = N ■ (58) 

Proof: See Appendix B. 
Theorem 3: For a MIMO broadcast system with N t transmit antennas and N r single- 
antenna receivers (N t > N r ) and fixed rates R\, Rn t , using the lattice-basis-reduction 
method, 

lim ~ lQgPe = N t . (59) 

p^oo log p 

Proof: Consider B = [bi...bjv r ] as the LLL-reduced basis for the lattice generated 
by H + . Each transmitted vector s is inside the parallelotope, generated by ribi, ....,r Nr h Nr 
(where r 1 ,...,r Nr are constant values determined by the rates of the users). Thus, every 
transmitted vector s can be written as 

s = ttbx + ... + t Nr b Nr , ^y<U<^. (60) 

For each of the transmitted vectors, the energy is 

P= ||s|| 2 = ||t 1 bi + ... + t JV - r bj Vr || 2 (61) 

^P<(\\tih 1 \\ + ... + \\t N h Nr \\) 2 (62) 

r lin n , , r N r 



P< (-^||bi|| + ... + -^||bjv r ||) . 



Thus, the average transmitted energy is 
Pav = E(P) < Ar^ max {^||b 1 ||,...,^||b JV J|}) 2 <c 1 .(max{||b 1 || 2 ,..,||b A r r || 2 }) (64) 

where c\ — ^f - max {r 2 , r% }. The received signals (without the effect of noise) are points 
from the l? Nr lattice. If we consider the normalized system (by scaling the signals such that 
the average transmitted energy becomes equal to one), 

d2 = PZ~ C 1 .(m aX {||b 1 ||2,...,||b JVr P}) (65) 
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is the squared distance between the received signal points. 

1 1 
For the normalized system, - is the energy of the noise at each receiver and — is the 

P 2p 
energy of the noise per each real dimension. Using (1651) . for any positive number 7, 



Using lemma 2, 



Pr I d 2 < 1 



< mTp iiTTlrn^ 

ci max{||bi|| 2 , ||bjv r |r} 



max{||b 1 ||, ||bjv-J|} < ^ (67) 



I =^ Pr <j d z < -L \ < Pr { — ^- < ±| = Pr < j (68) 



The JV( x jV r matrix H H has independent complex Gaussian distribution with zero mean 
and unit variance. Therefore, by using lemma 3 (considering M = N r and iV = N t ), we can 
bound the probability that <i H H is too small. 

Case 1, N r = N t : 



Pr jd 2 < 1 j < Pr |4„ < }> (69) 

2 \ N t r / 1 2 \ JVt+i 

7Cl«Ar t \ I / 1 7Cl«Ar t X 



(2 \ iV t 1 
1CiaNt ) maxlflnp)^ 4 " 1 ,!) for 7 > 1 and p > =- (71) 
p J I ) cio^ t 

< C2 ?[ (In p)^ -1 " 1 for 7 > 1 and p > max 1 e > (72) 



where c 2 is a constant number and e is the Euler number. 

If the magnitude of the noise component in each real dimension is less than ^d, the 
transmitted data will be decoded correctly. Thus, we can bound the probability of error by 
the probability that |u>j| 2 is greater than ^d 2 for at least one i, 1 < i < 2N t . Therefore, using 
the union bound, 
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P e <2N t [Prl l \w 1 \ 2 >-d< 



2N t ( Pr{d 2 < ^|.Pr|| Wl | 2 > Kl 2 



(73) 



d 2 < - 
P 



{-<d 2 <-\.Yi{\ Wl \ 2 >U 2 



P 



P 



4 l9 8 
- < d 2 < - 
P P 



8 l2 16 
- < d 2 < — 

P P 



(74) 



<2iV,(Pr^ 2 <l} + Pr{l<^<^.Pr{K| 2 >^ 
-<d?< —\ .Pr ( K| 2 > ^ .- 



< 2iV( ( Pr <jc/ 2 < -\ +Pr (d 2 < -l.Pr(| Wl | 2 > - 



(75) 



P 



P 



PHd 2 <^j.Pr jK| 2 >^ 



(76) 



For the product terms in ( 1761 ), we can bound the first part by (|72l . To bound the second 
part, we note that wi has real Gaussian distribution with variance Therefore, 



Pr<! M 2 > - p }=Q(V29)<e- 



(77) 



Now, for p > max ej, 

dH, (j76j) and (jTTD P e < 27V t (Vr || Wl | 2 > ^d 2 



2^(»+3) 



c 2 



N, 



(lnp) 



iV t +l _2* 



< 



(lnp) 



iVt+l 



,c 2 .2iV t 4 JVi + ^2^( J+3 ) e - 21 

\ i=0 

c 3 (\np) Nt+1 



< 



P 



(78) 
(79) 
(80) 
(81) 



where c 3 is a constant number which only depends on N t . Thus, 
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lim " lQgPe > lim ^loSP-Mlnp)^ +1 - logos = ^ 

According to Theorem 2, this limit can not be greater than N t . Therefore, 

lim Zi^ = Nt . (83) 

p^oo log p 

Case 2, N r < N t : 

For the N t x N r matrix H H , we use the first inequality in lemma 3 (by considering 
M = N r and N = N t ) to bound the probability that g? H h is too small: 

Pr|d 2 <2\< Pr (4 H <2£l<j (84 ) 



</W,(^f (85) 

< Pn, N r 1±Jk for 7 > 1 and p > =- (86) 

VP/ cia^ r 

C 2 l Nt 1 

< ^4— for 7 > 1 and p > =- (87) 

P Nt cia 2 Nr 

The rest of proof is similar to the case 1. 

■ 

Corollary 1: Perturbation technique achieves the maximum precoding diversity in fixed- 
rate MIMO broadcast systems. 

Proof: In the perturbation technique, for the transmission of each data vector u, among 
the set {H + (u + al)||l G Z 2Afr }, the nearest point to the origin is chosen. The transmitted 
vector in the lattice-reduction-based method belongs to that set. Therefore, the energy of the 
transmitted signal in the lattice-reduction-based method can not be less than the transmitted 
energy in the perturbation technique. Thus, the average transmitted energy for the perturbation 
method is at most equal to the average transmitted energy of the lattice-reduction-based 
method. The rest of the proof is the same as the proof of theorem 3. ■ 
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B. Fixed sum-rate 

When the sum-rate R sum is given, similar to the previous part, we can define the outage 
probability as the probability that the sum-capacity of the broadcast system is less than R sum . 

Theorem 4: For a MIMO broadcast system with N t transmit antennas, N r single- 
antenna receivers, and a given sum-rate R SU m, 

lim - X °Z P out < NtN ^ (g8) 

p->oo log p 

Proof: 

For any channel matrix H, we have [10] 

C sum = sup log | I Nr + pH H DH | (89) 

D 

where D is a diagonal matrix with non-negative elements and unit trace. Also, [29] 

, (2/>trH H DH) JVr (2ptrH H H) JVr 
|2pH H DH| < 1 N N r 1 = 1 N N r ] • ( 90 ) 

The entries of H have iid complex Gaussian distribution with unit variance. Thus 
trpH H H is equal to the square norm of an N t N r -dimensional complex Gaussian vector 
and has a chi square distribution with 2N t N r degrees of freedom. Thus, we have (similar to 
the equations I22H281 in the proof of theorem 2), 

p,^^,^ (91) 

Pr { trH H H < — ^ }> (92) 



* ptd (forp>^-^) (93) 

where c is a constant number. Now, 



l im ~^ p out = (94) 

p^oo log p 

Hm - logPr {sup D log |Ijy r + pH H DH| < R sum } 

p->oo log p 

< Hm - log Pr{sup D log |2pH H DH| < R sum } 

~ P->°o logp 
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By using ([9(1, (EH), and d%]): 



lim 



-logP c 



out 



< lim 



-log 



N t N r . 



(97) 



logp 



p— »oo 



logp 



The slope A^A^ for the SER curve can be easily achieved by sending to only the 
best user. Similar to the proof of theorem 3, the slope of the symbol-error rate curve is 
asymptoticly determined by the slope of the probability that | h max | is smaller than a constant 
number, where h max is the entry of H with maximum norm. Due to the iid complex Gaussian 
distribution of the entries of H, this probability decays with the same rate as p~ NtNr , for 
large p. However, although sending to only the best user achieves the optimum slope for the 
SER curve, it is not an efficient transmission technique because it reduces the capacity to 
the order of log p (instead of N r log p). 



Figure 2 presents the simulation results for the performance of the proposed schemes, 
the perturbation scheme [1], and the naive channel inversion approach. The number of the 
transmit antennas is N t — 4 and there are N r = 4 single- antenna users in the system. The 
overall transmission rate is 8 bits per channel use, where 2 bits are assigned to each user, 
i.e. a QPSK constellation is assigned to each user. 

By considering the slope of the curves in figure 2, we see that by using the proposed 
reduction-based schemes, we can achieve the maximum precoding diversity, with a low 
complexity. Also, as compared to the perturbation scheme, we have a negligible loss in the 
performance (about 0.2 dB). Moreover, compared to the approximated perturbation method 
[14], we have about 1.5 dB improvement by sending the bits, corresponding to the shift vector, 
at the beginning of the transmission. Without sending the shift vector, the performance of 
the proposed method is the same as that of the approximated perturbation method [14]. The 
modified perturbation method (with sending two shift bits for each user) has around 0.3 dB 
improvement compared to the perturbation method. 

Figure 3 compares the regularized proposed scheme with V-BLAST modifications of 
Zero-Forcing and Babai approximation for the same setting. As shown in the simulation 
results (and also in the simulation results in [1] and [14]), the modulo-MMSE-VBLAST 
scheme does not achieve a precoding diversity better than zero forcing (though it has a 
good performance in the low SNR region). However, combining the lattice-reduction-aided 
(LRA) scheme with MMSE-VBLAST precoding or other schemes such as regularization 
improves its performance by a finite coding gain (without changing the slope of the curve 



VII. Simulation Results 
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Fig. 2. Symbol Error Rate of the proposed schemes, the perturbation scheme [1], and the naive channel inversion approach 
for Nt = 4 transmit antennas and N r = 4 single-antenna receivers with the rate R = 2 bits per channel use per user. 



of symbol-error-rate). Combining both the regularization and the shift vector can result in 
better performance compared to other alternatives. 

Figure 4 compares the performances of the fixed-rate and the variable-rate transmission 
using lattice-basis reduction for N t = 2 transmit antennas and N r = 2 users. In both cases, 
the sum-rate is 8 bits per channel use (in the case of fixed individual rates, a 16QAM 
constellation is assigned to each user). We see that by eliminating the equal-rate constraint, 
we can considerably improve the performance (especially, for high rates). In fact, the diversity 
gains for the equal-rate and the unequal-rate methods are, respectively, N r and N t N r . 

VIII. Conclusion 

A simple scheme for communications in MIMO broadcast channels is introduced which 
is based on the lattice reduction technique and improves the performance of the channel 
inversion method. Lattice basis reduction helps us to reduce the average transmitted energy 
by modifying the region which includes the constellation points. Simulation results show that 
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Fig. 3. Comparison of the regularized proposed scheme with V-BLAST modifications of Zero-Forcing and Babai 
approximation (for Nt = 4 transmit antennas and N r = 4 single-antenna receivers with the rate R — 2 bits per channel 
use per user.). 



the performance of the proposed scheme is very close to the performance of the perturbation 
method. Also, it is shown that by using lattice-basis reduction, we achieve the maximum 
precoding diversity with a polynomial-time precoding complexity. 

Appendix A 

In this Appendix, we compute the second moment of a parallelotope whose centroid is 
the origin and its edges are equal to the basis vectors of the lattice. 

Assume that A is an M-dimensional parallelotope and X is its second moment. The 
second moment of \A is (|) X. The parallelotope A can be considered as the union 
of 2 M smaller parallelotopes which are constructed by ±|bi, ±|b 2 , ±|bj^, where b i5 
1 < i < M, is a basis vector. These parallelotopes are translated versions of \A with the 
translation vectors T, = ±|bi ± |b2 ± .... ± |bjvf, 1 < i < 2 M . The second moments of 
these parallelotopes are equal to (t^) M+ ~ X + ||Tj|| 2 Vol(|./l), 1 < i < 2 M . By the summation 



Sum-Rate = 8 bits per channel use, 2 transmit antennas, 2 users 



Symbol 

Error 

Rate 




SNR (dB) 



Fig. 4. Performance comparison between the fixed-rate and the variable-rate transmission for N t = 2 transmit antennas 
and N r = 2 single-antenna receivers with sum-rate 8 bits per channel use. 



over all these second moments, we can find the second moment of A. 



i=l 



2) x + WTtfwiqA) 



x+2 jw ^(iib 1 |i^+...+iib M in.voi(-^) 

^ + ^(l|bi|r + ... + ||b M || 2 ).Vol(^) 



X = -(||b 1 || 2 + ... + ||b A/ || 2 ).Vol(^l). 
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Appendix B 
Proof of Lemma 3 

Lemma 3 states that the probability that a lattice, generated by M independent N- 
dimensional complex Gaussian vectors, N > M, with a unit variance per each dimension, 
has a nonzero point inside a sphere (centered at origin and with the radius e) is bounded by 
Pn,m £2N for N > M, and (3 NM e 2N max {(- \ne) N+1 , 1} for N = M > 2. We can assume 
that e < 1 (for e > 1, lemma 3 is trivial because the probability is bounded). 



A. Case 1: M — I 

When M = 1, the lattice consists of the integer multiples of the basis vector v. If the 
norm of one of these vectors is less than e, then the norm of v is less than e. Consider 
the variance of the components of v as g 2 . The vector v has an iV-dimensional complex 
Gaussian distribution, / v (v). Therefore, the probability of this event is, 

-2JV 



Pr{||v|| < s} = f / v (v) rfv < f dv < (3. 

J||v||<£ ^]|v||<e " 6 

When the variance of the components of v is equal to one, we have, 



e 



(102) 



Pr{||v|| <e}<(3 NA e 21 \ (103) 



B, Case 2: N > M > 1 

Consider £( V;L ,...,v M ) as the lattice generated by V!,v 2 ,...,v A /. Each point of -^( Vl ,...,v M ) 
can be represented by V( zli ... jZM ) = Z\V\ + z 2 V2 + ••• + z M v M , where Z\, z M are complex 
integer numbers. The vectors vi,v 2 ,...,vjv/ are independent and jointly Gaussian. Therefore, 
for every integer vector z = ...,z M ), the entries of the vector V( 2lv . . i2M ) have complex 
Gaussian distributions with the variance 

qI=\\z\\ 2 q 2 = (\ Zl \ 2 + ... + \z M \ 2 )g 2 . (104) 
Therefore, according to the lemma for M = 1, 

27V 

Pr{||v( Zl ,.,„ 2M )||< £ }</3 Ar , 1 - (105) 

+ ... + \z M \ 2 ) 

Now, by using the union bound, 

Pr {d H <e}< X>{H v (*i,-,*m)II < e} (106) 

z^0 
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£ 2N 

z _^ U Z l| + ••• + \ Z M\ ) 



(107) 



0- 



NA 



E 

k l<||z[l<2 



_2/V 



I2JV 



.2N 



E 

2<||z||<3 



-2N 



\2N 



+ 



3<||z||<4 



I2JV 



(108) 



The M -dimensional complex integer points z = (zi, zm), such that fc < ||z|| < k + 1, 
can be considered as the centers of disjoint unit-volume cubes. All these cubes are inside the 
region between the 2M-dimensional spheres, with radii k— 1 and k+2. Therefore, the number 
of M-dimensional complex integer points z = (zi, ...,zm), such that k < ||z|| < k + 1, can 
be bounded by the volume of the region between these two 2M -dimensional spheres. Thus, 
this number is bounded by c 1 /c 2A/ ~ 1 for some constant ci. Therefore, 



E 



-2N 



1 27V 



< c\k 



2M-1 e 



2N 



fc<||z||<fc+l 



k 2N 



(109) 



(USED, (H09D Pr {d n <e}< Cl (3 N>l e 2N + 2 2M ~ 1 c 1 p_ 



NA 



-2N 



2 2N 



,2Af-l 



NA 



-2N 



3 2N 



(110) 



1 



fe=i 



J.2JV-2M+1 " 



(111) 



According to the assumption of this case, N > M; hence, 2N — 2M + 1 > 2. Therefore, 
the above summation is convergent: 



Pr{||v|| <e}<p N , M e 2N . 

'Throughout this proof, Ci, C2, ... are some constant numbers. 



(112) 



26 




radius — £ 
center — 



Fig. 5. The family of spheres «S V , 2 



C. Case 3: N = M > 1 

Each point of L( Vl ,..., Vjv ) can be represented by zvn — v, where v belongs to the lattice 
£(vi,...,vjv_i) and z is a complex integer. Consider S v as the sphere with radius e and centered 
at v. Now, zw N — v belongs to S iff the z\ N belongs to S v . Also, the sphere S v includes 
a point zvn iff <S V , Z includes v , where S v>z = — is the sphere centered at v/z with radius 

— (see figure 4). Therefore, the probability that a lattice point exists in S v is equal to the 

\z\ 

probability that v N is in at least one of the spheres {S VjZ }, z ^ 0. 

If we consider <i H as the minimum distance of L( vli ... iVJV ) and R as an arbitrary number 
greater than 1 : 

Pr {d H < e} = Pr { (L (V1) ... )V „) - 0) n 5 ^ 0} = Pr i ViV G |J |J <S V , 2 I (113) 
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<Pr<v*e |J Sv.J+PrWG |J S V> A (114) 

[ \m< R J I ii7ii> fl J 

In the second term of (|1 14J) . all the spheres have centers with norms greater than R and radii 
less than 1 (because \z\ > 1). Therefore, 

\J S v , z C {x| ||x|| > R- 1} (115) 
llfll>R 

dml) (HUD < Pr i v w G |J <S V ,* I +Pr{|| Vj v|| > fl-l} (116) 

I \m<R J 



<Pr<^v,ve |J <S v , z >+PrWe (J S v , z \ +Pr {|| Vjv || > R - 1} . 

(117) 

We bound the first term of (II 171) as the following: 

PrWe |J S v , z } (118) 

||f ||<fl,|*|<el-" 
V ||v||<2ii|«|>l 2R<\\v\\<3R \z\>2 

• ••+ ]T ^ Pr{ V7Y e 5 VjZ } j . (119) 



Noting that the pdf of \ N is less than or equal to \, 



fli<^ E ^voi( l s v , 2 )+ E Vo1 ^ 

yi|v||<2_R \z\>l 2R<\\v\\<3R \z\>2 



+ E Vol(5 Vj2 ) ) (120) 

Le 1 - Ar J J R<Hv||<(Le 1 - JV J+l) J R \z\>[e 1 - N \ 



1 / ^JV c 2iV 

- % N I Z^ Z^ \ Z \2N T Z^< Z^< l^pW 
||v||<2i? |z|>l ' ' 2_R<||v||<3iJ |z|>2 ' ' 
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E E 

Le 1 - JV Jil<||v||<(L£ 1 - JV J+l)fl |z|>le 1-JV J 

By using (11091 ), for one-dimensional complex vector z = z, 



c 2 e 



2N 



\2N 



(121) 



E 



\z\>i 



c 2 e 



2N 



I2JV 



5> E 

fc=i fc<|z|<fe+l 



-27V 



I2JV 



< 



C1C2S c^e 



E 

k=i 



fc2N-l — j2N-2 



(122) 



Now, 



una una < ^ ( £ 1 

\l|v||<2Ji 



E 

2_R<||v||<3 J R 



c 3 e 



2/V 



)2N-2 



+ 



C 3 £ 



2iV 



[e 1 - N iR<\\v\\<(le 1 - N i+l)R 1 J 



(123) 



Assume that the minimum distance of L( V i,...,v W -i) * s ^/v-i- The spheres with the radius 
djv-i/2 and centered by the points of -^( Vl ,...,v JV _ 1 ) are disjoint. Therefore, the number of 
points from the (N — l)-dimensional complex lattice L( V1 r .. )V ^_ 1 ), such that ||v|| < 2R, is 



bounded by 



c 4 (2R+d N _ 1 /2) 2 



(it is bounded by the ratio between the volumes of (2N — 2)- 



dimensional spheres with radii 2R + d^-i/2 and djv-i/2). Also, the number of points from 
^(vi,...^-!). such that (k-l)R < ||v|| < kR, is bounded by ^(^"^(ff+'fr-i) (it is bounded 



by the ratio between the volumes of the region defined by (k — 1)R — d^^i/2 < ||x|| < 
kR + (Im-i/2 and the sphere with radius d N _i/2): 



(H23]) < 



c 5 (2/? + d jV _ 1 /2) 2 jV - 2 

«AT_i 



-.5 



2iV 



I e l — ■ JVJ 

C5-R ~ (R + dN-l) 2N \^ 1 
d N-l k=2 K 



(TO < 



c 5 (2i? + (i J v_ 1 /2) 

>2/V-2 

«7V_1 



2N-2 



.e 2N + 



C5-R (i? + dx-i] 

WjV-l 



e 27V .ln( £ 



l-Ns 



< CqE 2N . max 



2N-2 



j2N-2 ' 
U 7V-1 



1 . max {— hie, 1} . 



(124) 



(125) 



(126) 



According to the proof of the case 2, we have Pr {d N ^i < 77} < I3n,n-iT] 2N ■ Therefore, 

(127) 



R 2N ~ 2 
E djv „ 1 <( max I , 9 ^_ 9 , 1 



«JV-1 
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< 1. Pr {djv-i > R} + 2 2N ~ 2 . Pr \ -R < d N _ x < R \+3 2N ~ 2 . Pr \-R < cZjv-i < -R 



2 

(128) 

< 1 + 2 2N ~ 2 . Pr {d N ^ <R} + 3 2N ~ 2 . Pr jtfo-i < ii? J + ... (129) 

fc=i 

Erfjv_i ^ cq£ 2JS . max ( £.n-2 > ^ ) ■ max { — m£ > 1} r — °8£ 2N -R 2N ■ max { — lne, 1} . 

(131) 

To bound the second term of (II 171) . we note that for |z| > e 1 "^, the radii of the spheres 
S v z are less or equal to e N , and the centers of these spheres lie on the (N — 1) -dimensional 
complex subspace containing L vli ... iVjv _ 1 . Also, the norm of these centers are less than R. 
Therefore, all of these spheres are inside the region A which is an orthotope centered at the 
origin, with 2N real dimensions (see figure 5): 

U S v ,z C A (132) 

\\^\\<R,\z\>e 1 - N 

=}► PrW e |J S V J < -^Vol I |J S v , z < ^Vol(A) (133) 

< \{2e N ) 2 (2R + e N ) 2N - 2 (134) 

Also, according to the Gaussian distribution of the entries of \ N (which have Variance 
\ on each real dimension), we can bound the third term of (11171) as, 



By using (TOTT) . (fT34l) . and (fT35l) . 



Pr{d H <e} < c 8 £ 2Ar . J R 2Ar .max{-ln£,l} + 47(2e Ar ) 2 (2i?+£ Ar ) 27V - 2 +c 9 e~(^ . (136) 
The above equation is true for every R > 1. Therefore, using i? = v^^a/ - hi(e 2Ar ) + 1, 

Pr{rf H < e} < /3Ar,jv£ 2iV .max{(-ln£) JV+1 ,l} . (137) 
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Sul ispacc containing L, 



Fig. 6. The orthotope .4 
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